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The 3’ end of the 20-kb genome of the Purdue strain of porcine transmissible gastroen- 
teritis coronavirus (TGEV) was copied into cDNA after priming with oligo(dT) and the 
double-stranded product was cloned into the PsiI site of the pUC9 vector. One clone of 
2.0-kb contained part of the poly(A) tail and was sequenced in its entirety using the 
chemical method of Maxam and Gilbert. Another clone of 0.7 kb also contained part of 
the poly(A) tail and was sequenced in part to confirm the primary structure of the most 
3’ end of the genome. Two potential, nonoverlapping genes were identified within the 3'- 
terminal 1663-base sequence from an examination of open reading frames. The first gene 
encodes a 382-amino acid protein of 43,426 mol wt, that is the apparent nucleocapsid 
protein on the basis of size, chemical properties, and amino acid sequence homology with 
other coronavirus nucleocapsid proteins. It is flanked on its 5’ side by at least part of the 
matrix protein gene. The second encodes a hypothetical 78-amino acid protein of 9101 mol 
wt that is hydrophobic at both ends. A 3'-proximal noncoding sequence of 276 bases was 
also determined and a conserved stretch of 9 nucleotides near the poly(A) tail was found 
to be common among TGEV, the mouse hepatitis coronavirus, and the avian infectious 


bronchitis coronavirus. 


INTRODUCTION 


The genome of the porcine transmissible 
gastroenteritis coronavirus (TGEV) has 
been shown to be a single-stranded, non- 
segmented, polyadenylated, infectious RNA 
molecule of 6.8 < 10° mol wt or approxi- 
mately 20 kb in length (Brian et al., 1980). 
The total number of genes encoded by the 
TGEV genome, however, has not yet been 
determined. The genome codes for at least 
four unique polypeptides on the basis of 
existing protein data. The virion is com- 
prised of three major structural proteins: 
a 200-kd peplomeric glycoprotein, a 29-kd 
membrane-associated matrix glycoprotein, 
and an internal phosphorylated nucleo- 
capsid protein that measures from 46 to 50 
kd (Garwes and Pocock, 1975; Moreau and 
Brian, unpublished). These proteins alone 
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would account for only approximately 8.4 
kb of coding information. In addition, the 
virus synthesizes at least one nonstructural 
protein during its replication, an RNA-de- 
pendent RNA polymerase, the size of which 
is not yet known (Dennis and Brian, 1982). 
During replication, TGEV produces nine 
species of subgenome-size polyadenylated 
RNA moiecules each of which may function 
as a separate mRNA (Dennis and Brian, 
1982), assuming that the 3’ coterminal 
“nested set” arrangement described for the 
mRNAs of mouse hepatitis virus (MHV; 
Lai, et al, 1981; Leibowitz et al., 1982; Rot- 
tier et al, 1981) and the avian infectious 
bronchitis coronavirus (IBV; Stern and 
Kennedy, 1980) is also true for TGEV. 
From this information, TGEV may code for 
as many as 10 different protein species. 
One powerful approach for determining 
the number of potential genes in an RNA 
virus genome is to examine the primary 
nucleotide structure and deduce the iden- 
tity of genes from an examination of open 


reading frames. In this paper we describe 
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experiments that begin to examine the 
TGEV genome by this approach. cDNA 
clones were prepared from the 3’ terminal 
10% of the polyadenylated genome and 
were sequenced. Two potential genes were 
identified within the first 1663 bases. One 
gene encodes a protein of 382 amino acids 
which is the apparent nucleocapsid protein 
on the basis of size, chemical properties, 
and significant amino acid sequence ho- 
mology with other coronavirus nucleocap- 
sid proteins. This gene is flanked on its im- 
mediate 5’ side by at least part of the ma- 
trix protein gene. The second gene lies to 
the 3’ side of the nucleocapsid protein gene 
and encodes a hypothetical protein of 78 
amino acids that is hydrophobic at both 
ends. A 3’ noncoding sequence of 276 bases 
sharing a 9-base conserved sequence near 
the poly A tail with other coronaviruses 
was also determined. 


MATERIALS AND METHODS 


Virus and cells. The Purdue strain of 
TGEV was plaque purified and grown on 
the swine testicle (ST) cell line as previ- 
ously described (Brian et al., 1980). 

Purification of genomic RNA. Virus was 
purified from clarified supernatant fluids 
as previously described (Brian et al., 1980) 
except that all sucrose solutions were made 
up in TMEN (10 mM Tris-maleate, pH 6.0, 
100 mM NaCl, 1 mM EDTA). Viral RNA in 
1 of 10 flasks was radiolabeled in order to 
follow the purification of the RNA. For 
these radiolabeled cultures, infected cells 
were refed with phosphate-free medium 
containing 1% fetal calf serum and 40 uCi 
[*P]orthophosphate (ICN) per milliliter. 
Viral RNA was extracted by dissolving the 
virus pellet in 0.5 ml TNE (10 mM Tris- 
hydrochloride, pH 7.5, 100 mM NaCl, 1 mM 
EDTA) containing 1% SDS and 0.5 mg pro- 
teinase K per milliliter, incubating for 0.5 
hr at 37°, and extracting twice with an 
equal volume of a mixture of 50% phenol/ 
48% chloroform/2% isoamyl alcohol. RNA 
was ethanol precipitated after adding 0.1 
volume 2 M sodium acetate. Because small 
molecular weight RNA species are found 
in some preparations of purified corona- 
virion RNA, full-length genomic RNA to 


be used for cDNA cloning and making 
probe for colony screening was selected by 
rate-zonal sedimentation on preformed 
linear gradients of 15 to 30% sucrose (wt/ 
wt) made up in TNE-0.1% SDS. RNA was 
dissolved in water and sedimented 1.5 hr 
at 110,000 g, 25°, on 5-ml gradients. Frac- 
tions of 0.2 ml were collected and the dis- 
tribution of radioactivity was determined 
by Cerenkov counting. Only RNA sedi- 
menting with a sedimentation coefficient 
of 50 S or greater, as determined by ref- 
erence to sedimentation of mammalian 28 
S and 18 S ribosomal RNA in a parallel 
gradient, was recovered by ethanol precip- 
itation and used in the experiments de- 
scribed below. 

cDNA cloning of the 3’ end of the TGEV 
genome. TGEV genomic RNA was cloned 
using a modified method of Gubler and 
Hoffman (1983). First strand synthesis was 
carried out in a reaction volume of 50 ul 
containing 50 mM Tris-hydrochloride, pH 
8.3, 10 mM MgCl,, 10 mM DTT, 2 mM dCTP, 
2 mM dTTP, 2 mM dATP, 2 mM dCTP, 10 
uCi [P]dCTP (3000 Ci/mmol, ICN), 50 
pmol oligo dTj2-13, 6 ug TGEV RNA, 30 U 
RNAsin, 10 U reverse transcriptase (Sei- 
kagaku), for 1 hr at 42°, and the reaction 
was stopped by adding 2 wl 0.5 M EDTA. 
Nucleic acids were phenol-chloroform- 
isoamy] alcohol extracted and ethanol pre- 
cipitated after the addition of 0.5 vol of 7.5 
M ammonium acetate. 

Second strand synthesis was carried out 
in a reaction volume of 100 yl containing 
20 mM Tris-hydrochloride, pH 7.5, 5 mi 
MgCle, 10 mM (NH,)2S0,, 100 mM KCl, 0.15 
mM B-NAD, 50 mg/ml BSA, 40 nM dNTPs, 
8.5 U/ml Escherichia coli RNAse H, 230 U/ 
ml DNA polymerase I, 10 U/ml DNA li- 
gase, and all of the product from the first 
strand reaction. The reaction was incu- 
bated at 12° for 1 hr, then at 22° for 1 hr. 
The reaction was stopped by adding 4 ul 
0.5 M EDTA and reaction products were 
phenol-chloroform-isoamyl alcohol ex- 
tracted and fractionated on a Sephadex 
G50 spun column (Maniatis et al, 1982), and 
the ds cDNA was ethanol precipitated. 

Double-stranded cDNA was homopoly- 
mer tailed essentially by the method of 
Roychoudhury and Wu (1980). The follow- 
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ing were added together, in order: 3 pl 
daCTP” (>800 Ci/mmol), 20 ul 10X caco- 
dylate buffer (1.4 mM potassium cacodyl- 
ate, 0.3 M Tris-hydrochloride, pH 7.6), 4 ul 
5 mM DTT, 3 «1 10 mM dCTP, 2 ul 100 mM 
CoCl,, 12 units terminal deoxynucleotide 
transferase (PL Biochemicals; at least 8 
units/pmol 3’ end) in 1.5 ul, H,O to 200 ul 
final volume. The reaction was carried out 
at 12° for 1.5 min then stopped by adding 
10 wl 0.5 M EDTA. This reaction resulted 
in an average of 15 dCMP residues added 
per 3’ end of dsDNA, the optimal number 
for annealing and transformation (Peacock 
et al., 1981). 

C-tailed ds cDNA was annealed to G- 
tailed, PstI-linearized pUC9 vector (PL 
Biochemicals) for 4 hr at 58° in a 50-1 vol 
of buffer containing 10 mM Tris-hydro- 
chloride, pH 7.5, 150 mM NaCl, 1 mM 
EDTA. The total concentration of DNA 
used was less than 0.5 pg/ml and the op- 
timal insert:vector ratio was 1:1 on a mass 
basis. F. coli strain JM103 was transformed 
using the method of Hanahan (Hanahan, 
1983). Cells containing inserts were ob- 
served as white colonies on YT agar plates 
that contained 100 pg ampicillin/ml, 1 mM 
IPTG, and 0.004% X-gal (Heidecker and 
Messing, 1983). Recombinant colonies were 
transferred to nitrocellulose (Millipore, 
HAWP) and probed with random-primed 
cDNA copied from TGEV genomic RNA. 

Identification of large clones containing 
3'-specific TGEV sequences. “P-labeled, 
random-primed cDNA used for colony hy- 
bridization was synthesized as described 
above for the oligo(dT)-primed reaction 
except that 0.2 ug of RNA was used and 
oligo(dT) was replaced by 20 ug of frag- 
mented calf thymus DNA. Probe was alkali 
treated to hydrolyze the RNA and then was 
used for colony hybridization (Maniatis et 
al., 1982). Colonies yielding a strong sig- 
nal were analyzed for insert size by elec- 
trophoresis of plasmid DNA in agarose 
gels (Kado and Liu, 1981). Inserts of 0.2 to 
2.0 kb (the largest) were further analyzed 
by Southern hybridization with *P-la- 
beled poly(dT) to detect poly(dA) content 
and by cross-hybridization with nick- 
translated inserts to detect overlapping 
sequences. **P-labeled poly(dT) probe was 


prepared as described above for the 
oligo(dT)-primed reaction except that 50 
pmol oligo(dT): poly(rA) (PL Biochemi- 
cals) replaced the RNA. Alkali-treated ”P- 
poly(dT) probe was incubated for hybrid- 
ization at 37° for 12 hr then at 20° for 36 
hr, and blots were washed in 2X SSC, 0.1% 
SDS at 20°. 

Restriction endonuclease mapping. Plas- 
mid was purified by lysozyme lysis and ce- 
sium chloride centrifugation (Maniatis et 
al, 1982), and restriction endonuclease 
mapping was done essentially as described 
by Smith and Bernstiel (1976) using plas- 
mids that were labeled at the SalI site 
within the multiple cloning linker region. 

DNA sequencing and sequence analysis. 
Restriction fragments end labeled with “P 
were isolated and sequenced by the method 
of Maxam and Gilbert (1980). Sequences 
were analyzed with the aid of the program 
developed by Queen and Korn (1984) and 
sequence homologies were searched against 
Genbank, both marketed as part of the 
Beckman Microgenie program, March 1985 
version (Beckman Instruments, Inc.). 


RESULTS 


cDNA cloning and sequencing of two 
clones from the 3’ end of the genome. Start- 
ing material for cDNA cloning was ap- 
proximately 6 wg of rate-zonally purified 
genomic RNA obtained from 400 ml of tis- 
sue culture medium. An estimated 200 ng 
of ds cDNA was obtained, as determined 
by radiolabel incorporation during second 
strand synthesis, and from this approxi- 
mately 2000 white colonies were obtained. 
By colony screening 200 colonies gave a 
strong signal to “P-labeled cDNA prepared 
from genomic RNA, and of these, 13 had 
inserts of 200 to 2000 bases as determined 
by agarose gel electrophoresis of super- 
coiled plasmids, and were further analyzed 
by restriction enzyme analysis and poly (A) 
content. The largest clone of 2000 bases, 
FG5, did not react by Southern blotting to 
P-labeled oligo(dT), but did cross-hy- 
bridize in Southern blot analysis with sev- 
eral other smaller clones that did react 
strongly with oligo(dT). One of these, J21, 
a clone of 700 bases, was sequenced in part 
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to determine the primary structure of the 
extreme 3’ end of the genome. 

The orientation of clones FG5 and J21 in 
reference to the virus genome and the re- 
striction enzyme sites used for sequencing 
are illustrated in Fig. 1. Our orientation 
presumes polyadenylation at only the 3’ 
end of the genome and this, in turn, is based 
on the precedent of the documented 3’ 
polyadenylation site in the avian infectious 
bronchitis virus and mouse hepatitis virus 
genomes (Lai et al, 1981; Stern and Ken- 
nedy, 1980). The strategy used for sequenc- 
ing is described in the legend to Fig. 1. Over 
96% of the sequence containing the two 
complete genes we report was determined 
either by sequencing both strands or by re- 
peated sequencing of the same strand using 
different methods of end labeling. Some of 
the sequences were derived from subclones 
of FG5 made from the PstI restriction sites. 

The total sequence of FG5 is illustrated 
in Fig. 2. Sequences from J21 thal overlap 
with FG5 are identical to those of FG5 ex- 
cept that the total length of the polyade- 
nylate tail is 15 bases for the J21 clone, and 
6 for the FG5 clone. 

The entire nucleotide sequence was 
translated in all possible reading frames 
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and only translation of the virus-sense 
strand revealed open reading frames of 
greater than 120 bases that are preceded 
by a termination codon and contain an ap- 
propriate initiator methionine codon (Fig. 
3). Of these, only the two largest open 
reading frames are evaluated below. 

The largest open reading frame predicts 
a protein having properties expected of the 
nucleocapsid protein. The largest open 
reading frame extends from base 353 to 
base 1498 and predicts a 382-amino acid 
protein of 43,426 mol wt. The only TGEV 
structural or nonstructural protein de- 
scribed to date that approaches this size is 
the phosphorylated nucleocapsid protein 
that measures 46 to 50 kd by SDS-poly- 
acrylamide electrophoresis (Garwes and 
Pocock, 1975; Moreau and Brian, unpub- 
lished). The protein has two properties that 
are strikingly similar to the nucleocapsid 
proteins of MHV and IBV (Armstrong et 
al., 1983; Boursnell et al., 1985; Skinner and 
Siddell, 1988). First, it is rich in serine. 
Thirty-nine (10% )of the residues are serine 
making it the most abundant amino acid. 
Assuming this protein is phosphorylated 
at serine residues, as is the MHV A59 pro- 
tein (Stohlman and Lai, 1979), then a high 


kb 
6 5 4 3 2 1 Og enq 


—— pTGEV J21 
pTGEV FGS 


= ae - ae 
e is £ a3 
oO 500 1000 1500 2000 
o—_—_—_—_—-- eo —____——_ OO —_—_——_ ————_—-e 
=——_——_ , e—_—_—_——__ o—_——_ 
Oo 
Se = ee 
_e ———_--e —_———-——- 


Fic. 1. Restriction endonuclease map and sequencing strategy for TGEV cDNA clones J21 and 
¥FG5. The internal HindIII, TaqI, Pst, and Xbal sites, derived by restriction endonuclease mapping, 
and the HindIII and Sail sites in the multiple cloning region of the pUC9 vector, were the sites used 
for initial DNA sequencing. Internal Aécl, Kpnl, and Ndel sites were identified from sequence data 
and were used to complete the sequencing. @ Indicates sites labeled at the 5’ end using polynucleotide 
kinase. @ Indicates sites labeled at the 3’ end using reverse transcriptase and the appropriately 
labeled deoxynucleotide triphosphate. @ Indicates site labeled at the 3’ end using dideoxy A and 


terminal transferase. 
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98 128 
cca CEST TCCOTTR ICE PTNGGAAGRAGCTATNCCTIOCTCTOGAASGTOTIOCAACT GGG ICACTCTAACTTOCTTTCAGCARPTIGTACECTGAAGG 
LCOCVSALGRSY¥YVLPLEGVPTGVTLTLELESGNLYAEG 


158 180 210 . . 246 
GTTCAAAATIGCAGCTGGTATGAACATCGACAATTTACCAAAATACGTAA TGGTTGCA TTACCTAGCAGGACTA TIGICTACACACTTGTIGGCAAGAAGTTGAAAGCAAGTAGTGCGAC 
FKIAGGMNIODBNLPKYVMVALPSRTIVYTEVGKKLKAS SAT 


. 270 . 348 . 338 360 
GCATGOSCTTACTATGTAAAATCTAAAGCTGGIGATT. (TTACTCAACAGAGGCAAGAACTGA TAA TPIGAGTGAGCAACAAAAATTATTACA TATGGTATAACTAAACTICTAAATGOOCAA 
GwaAyYYVKS kKAGOYSTEARTOWNWLS EQEKLUH MV MAN 
. . 396 . 420 456 480 


CCAGGGACAACCTGTCAGTTGGGGAGATGAATCTIACCA AAACACGT GG TOGTTOCAA PICOOGTGGTOGGAAGAATAATAACATACCTCTTTCATICTTCAACOOCATAACOCTOCAACA 
eceRVSWGCDESTKTRGRSNSRGRKNNNI PLS FFNPITL QQ 


51¢ 546 576 . 680 
AGGTTCAAAATTPTGGAACTTATGTOOGAGAGACT TTGTACOCAAAGGAATAGGTAACAGGGATCAACAGATIGGT'TATIGGAATAGACAAACTOGCTATOGCATGGTGANGGGOCAACG 
cs K FWNLCPRDEFVPKGIGNRBOQOQIGY WNROT RY RMVK GOR 


. . 636 668 698 728 
ITAAAGAGCPICCTGAAAGGTGGTICTICTACTACTTAGGTACTOGACCTCATIGCAGATGCCAAA TTTAAAGATAAATTAGA TGGAGTIGICTGGGTIGOCAAGGATGGTGOCATGAACAA 
K EL P ERWFFYY¥YLGTGPHAODAK FKODKLODGVVWVAK OGAMN K 


. 758 788 818 848 
AOCAACCAQACTIGGTAGTCGTGGTGCTAATAATGAATOCAAAGCTTTIGAAA TIOGA TIGGTAAAGTGOCAGGUGAA TITCAACTTGAAGT TAA TCAATCAAGAGACAATICAAGGTCACG 
PTT LGS R GANNES KALK FOGKVPGEFQLEVNQS RDN S RSER 


878 968 930 . . 968 
CTCICAATCTAGATCTOGGTCTAGAAATAGATCTCAATCTAGAGGCAGGCAACAA PICA TAACAAGAAGGATGACAGTGTAGAACAAGCIGTICTTIGOCCACTTAAAAAGTTAGGTGT 
sQ@QSRSRSRNRSQSRGRQQFNNKKODODSVEQAVLAALK KLUG VY 


. . 998 . 1628 . 1858 . . 1489 
‘TGACACAGAAAAACAACAGCAACGCTCTOGTICTAAATCTAAAGAACGTAGTAACTCT AAGACAAGAGATACTACACCTAAGAA TGAAAACAAACACACCTOGAAGAGAACTGCAGGTIAA 
DT EKQQQRSRSK SKERSNSKTROTTPK NEN KHTS KRTFAG K 


lla - - 1146 = - 1176 ° . 1296 
AGGTGATGTGACAAGATITTATGGAGCTAGAAGCAGTICAGOCAA TITTGGTGACACTGACCTOGTTIGOCAA TGGGACCAGTGCCAAGCA TTACCCACAACTGGCIGAATGTGTICCATC 
GDVTRFYGARSSSANFGOTODLVANGSSAKH Y P QLA ECV PS 


VS SILFGSsSYWTS KEDGOeQt1tEVTF TH K ¥Y HLPK DODPKTGOQFL 


. 135¢ . . 1388 > 141¢ . 1448 
ITCAGCAGA TAA TGOCTATGCICGTOCA TCAGAAGIIGGCAAAAGAACA GA GAAARAGAAAATCTOGPICTAAATCTGCAGAAAGGTCAGAGCAAGATGTGGTACCTGATGCATTAATAGA 
QQrNAYARPSEVAKEQRKRKSRSKSAERSEQODVVPODALIE 


. 1478 1586 . 1530 . 1560 
AARTTATACAGATGTGTTTGATGACACACAGGTTGAGATARTTGATGA GGTAACGAACTAAACCAGATOCTOGTCTTIOCIOCATGCIGTATTTATTACAGTTTTAATCTTACTACTAATT 

NYTODVFODODTQVETIODEVTN MLVFUHAVFITVLILLUI 
. 1596 . 1626 1658 . 1680 


GGTAGACTOCAATTATTAGAAAGACTATTACTIGATCACICTTTCAATCIT “AAAACTGICAATGACTTTAATATCITA TATAGGAGTTTAGCAGAAACCAGATTACTAAAAGTOGTGCTT 
GrRUbLUEQutuLERUtCLtCoOoH SFNOBKT VNODF NIT LY RS LAETREL KV VE 


RVvVIFLVLELGFEFCCYRELVTLV 


. 1838 . 1868 . 1898 . 1928 
ACTCTIGTACAGAATGGTAAGCACGTGTAA TAGGAGGTACAAGCAROOCTATIGCATATTAGGAAGTTTAGATTTGATTIGGCAATGCTAGA TTTAGTAATTTAGAGAAGTTTAAAGATC 


195¢ . 1988 . 2018 
CGCTACGACGAGRCAACAATGGANGAGCTAACGTCTOGATCTAGTGATTGTTTAAAATGTAAAATIGITIGAAAATTTTOCTTTTGATAGTGA TACAAAARACOOOOCOOOOOKL 


Fic. 2. The primary nucleotide sequence of clone FG5 and the deduced amino acid sequences for 
a portion of the matrix glycoprotein (bases 20 through 387 in the second reading frame), the nu- 
cleocapsid protein (bases 353 through 1498 in the second reading frame), and the hypothetical hy- 
drophobic protein (bases 1507 through 1740 in the first reading frame). A 10-base sequence highly 
conserved among coronaviruses is underlined beginning at base 1940. 


level of phosphorylation might explain the only 46 (12% )are acidic, giving the protein 
3- to 6-kd difference between the predicted a net charge of +28 at neutral pH. 

and measured molecular weights. Second, Although the consensus sequence around 
the protein is basic, a property expected of the AUG initiator codon for the TGEV nu- 
nucleic acid-binding proteins. Sixty-nine cleocapsid protein (UAAAUGG) is not 
(18% )of the amino acids are basic whereas among the most favored for translation 
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Fic. 3. Schematic diagram of possible open reading frames obtained when translating the FG5 
nucleotide sequence as either virus-sense RNA or virus complementary-sense RNA. Vertical bars 
above the line represent the first methionine codon that could serve as the initiation site for trans- 
lation. Vertical bars below the line represent termination codons. (M), partial sequence of the matrix 
protein gene. N, sequence of the nucleocapsid protein gene. HP, sequence of a hypothetical protein 


gene. 


initiation sites that have been described, 
it is not without precedent (Kozak, 1983). 
This AUG, therefore, probably identifies 
the authentic beginning of the TGEV nu- 
cleocapsid gene since the sequence from 
this point leftward to the end of clone FG5, 
except for a 12-base intergenic sequence, 
reveals an open reading frame coding for 
a protein sharing extensive regions of 
amino acid homology with the small matrix 
glycoprotein (M or E1) of the mouse hep- 
atitis virus A59 (Armstrong et al, 1984; 
discussed below). 

A second open reading frame to the 3’ side 
of the nucleocapsid protein gene encodes a 
hypothetical protein of 9101 mol wt that is 
hydrophobic at both ends. An open reading 
frame beginning at base 1507 and extend- 
ing through base 1740 encodes a hypothet- 
ical 78 amino acid protein of 9101 mol wt 
(Fig. 2). A hydrophobicity analysis of the 
protein reveals that it is hydrophobic for 
a stretch of approximately 25 amino acids 
at each end and it is hydrophilic in its cen- 
tral region. There are eight basic amino 
acids and four acidic amino acids giving 
the protein a net +4 charge at neutral pH. 
Basic and acidic amino acids are distrib- 
uted evenly throughout the central hydro- 
philic region, but 4 basic amino acids and 
no acidic ones are found among the 27 


amino acids at the carboxy terminus. There 
is yet no direct evidence for this protein. 


DISCUSSION 


We present the primary nucleotide se- 
quence for the TGEV nucleocapsid protein 
(N) gene and the deduced amino acid se- 
quence for the protein. This is the first pri- 
mary sequence data for a coronavirus in 
the antigenic subgroup to which TGEV be- 
longs, and such information allows one to 
first, firmly conclude that TGEV shares an 
ancestral relationship with MHV and IBV, 
and second, to identify potentially func- 
tional domains on the N protein by ex- 
amining conserved structures among the 
diverged viruses. The first two coronavirus 
N gene sequences to be described are those 
of the closely related JHM and A59 strains 
of MHV (Armstrong e¢ al., 1983; Skinner 
and Siddell, 1983) and between these an 
overall homology of 94% was found for both 
the nucleotide and amino acid sequences, 
reflecting the antigenic similarities be- 
tween the viruses. Interestingly, the anti- 
genically distinct avian infectious bron- 
chitis virus shows no N gene nucleotide se- 
quence homology with MHV, yet shares an 
overall amino acid sequence homology of 
26% (Boursnell et al., 1985). Furthermore, 
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there is a stretch of 67 amino acids within 
the amino terminal one-third of the protein 
that shows a sequence homology of 51% 
between the viruses (Boursnell et al, 1985) 
suggesting that a strong selective pressure 
exists for a specific functional group defined 
by this sequence. This interesting pattern 
repeats itself in the structure of the TGEV 
N protein. Although TGEV shows no an- 
tigenic relatedness to either MHV or IBV 
(Pedersen et al., 1978), and no N gene nu- 
cleotide homology with either MHV or IBV, 
it does show an overall amino acid homol- 
ogy of 27 and 26% with MHV (JHM) and 
IBV, respectively. Furthermore, the con- 
served 67 amino acid region is also found 
in TGEV (becoming 68 positions when 
TGEV is compared; Fig. 4). This conserved 
region is slightly more basic than the over- 
all nucleocapsid protein and therefore may 
function as a site of interaction with ge- 
nomic RNA. 

Other regions in the N proteins of the 
three viruses share structural similarities 
in the absence of a common primary struc- 
ture suggesting the existence of additional 
conserved functional domains. Although 
the N proteins are different lengths (382 
amino acids for TGEV, 455 for MHV, and 
409 for IBV), when the three are aligned 
by the 68-amino acid conserved sequence, 
the following structural similarities are 
observed. (i) Four cluster groups contain- 
ing 2-10 serine residues are found in par- 
allel with TGEV amino acid positions 20- 
40, 150-190, 260-300, and 340-360. Other 
smaller serine clusters are found in MHV. 
In all three viruses, regions of 10-40 amino 
acid stretches can be found that are void 
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Fic. 4. Amino acid sequence homologies among IBV, 
MHV, JHM, and TGEV for a 68-amino acid conserved 
region in the nucleocapsid protein gene. The sequence 
starts at amino acid 53 from the initiator codon for 
IBV (Boursnell et al., 1985), amino acid 86 for MHV 
JHM (Skinner and Siddell, 1983) and amino acid 53 
for TGEV. Identical amino acids are boxed in. 


of serines. (ii) Three cluster groups of 5- 
29 basic amino acid residues are found in 
parallel with TGEV amino acid positions 
0-30, 150-260, and 330-350. (iii) A cluster 
of 9-11 acidic amino acid residues is found 
within the last 32 amino acids at the car- 
boxy terminus. 

Although TGEV would appear to be 
equally diverged from IBV and MHV on 
the basis of amino acid sequence, TGEV 
more closely resembles MHV in its genome 
arrangement. Firstly, like MHV, the N 
gene for TGEV is flanked on its 5’ side by 
the matrix protein (M or E1) gene, whereas 
for IBV, two genes, derived from overlap- 
ping reading frames and encoding hypo- 
thetical proteins of unknown function, lie 
between the M and N genes (Armstrong et 
al., 1984; Boursnell and Brown, 1984). Our 
conclusion that the M gene for TGEV lies 
to the immediate 5’ side of the nucleocapsid 
gene is based on amino acid sequence ho- 
mology with the M protein of MHV A59. 
Of the 105 amino acids deduced for the 
TGEV matrix protein sequence, 31% are 
perfectly homologous and another 15% are 
conservative differences (Fig. 2 and data 
not shown). Secondly, the number of nu- 
cleotides separating the M (E1) and N 
genes is close, 14 for MHV and 12 for 
TGEV, and these match perfectly for a 
stretch of 8 bases: 


MHV A59 TCTAAACTTTAAGG 
TGEV CTAAACTTCTAA 


Since part of this sequence may play a role 
in primer recognition for transcription 
(Brown and Boursnell, 1984; Budzilowicz et 
al., 1985), some common features between 
the leader molecules of MHV and TGEV 
may be anticipated. 

No direct evidence exists for the hydro- 
phobic protein encoded by base 1507 
through 1740. Genes encoding small hy- 
drophobic proteins in MHV and IBV have 
been described, however (Boursnell and 
Brown, 1984; Skinner ef al., 1985), but their 
hydrophobicity is only at one end, they map 
at an entirely different region in the ge- 
nome, and no sequence homology is found 
between them and the TGEV hydrophobic 
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protein. Regarding this open reading 
frame, it is noteworthy that three small 
polyadenylated, putative messenger RNAs 
have been identified in TGEV-infected cells 
that have not been reported for MHV or 
IBV (Dennis and Brian, 1982). Assuming 
TGEV replicates by the consensus scheme 
proposed for the replication of MHV and 
IBV, namely that all messages have a 3’- 
coterminal nested set arrangement (as 
suggested by preliminary experiments with 
TGEV [Hu et al, 1984]), then one of the 
small messages described by Dennis and 
Brian may be the message for the hydro- 
phobic protein. From the known sequence 
(Fig. 2), such a message would be 0.20 Md. 
Two structural features favor the plaus- 
ability of this being a functional hydro- 
phobic protein gene. (i) The intergenic se- 
quence preceding the gene, inclusive of the 
N gene stop codon, contains a 6-base se- 
quence, CTAAAC, that is in common with 
part of the intergenic sequence preceding 
the N gene for both MHV and TGEV de- 
scribed above, and may play a role in the 
initiation of mRNA transcription (Budzi- 
lowicz et al., 1985). (ii) The 7-base sequence, 
GAGAUGC, at the initiation site of the 
hydrophobic protein is a favored pattern 
among eukaryotic initiation sequences 
(Kozak, 1983). 

Assuming that the gene for the hydro- 
phobic protein is real then the 3’ terminal 
noncoding sequence would be a total of 276 
bases, exclusive of the poly (A) tail, and 
would be the shortest noncoding sequence 
of those identified for coronaviruses. The 
significance of the noncoding region is not 
completely known although it undoubtedly 
functions as an attachment region for the 
polymerase to initiate synthesis of the 
negative strand RNA. One possible site 
that may be critical for recognition or 
binding is a 10-base sequence, GGGAA- 
GAGCT, that is conserved between IBV 
(found 81 bases from the 3’ end) and MHV 
(82 bases from the end). With the exception 
of the first base, a T instead of G, TGEV 
shares an identical sequence beginning 77 
bases from the 3’ end (Fig. 2). 
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